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We demonstrate that a broad spectrum of spreading processes taking place on complex networks 
can be categorized into two fundamentally different regimes, in which either positive or negative 
degree-degree correlation decelerates their propagation. This result can be explained by the role 
of the nodes with a high connectivity and their distribution in the network, while they act either 
as accelerators or delayers. This finding is relevant for controlling the velocity of many different 
spreading processes, ranging from epidemic diseases to dissemination of information and cascading 
failures in power grids. 
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Understanding the mechanisms of spreading phenom- 
ena is a need shared across many scientific disciplines, 
with examples as seemingly diverse as reaction diffu- 
sion processes, pandemics or cascading failures in electric 
power grids. Substantial new insights have recently been 
gained through the application of statistical physics on 
the study of large-scale networked systems, the results 
of which show how the underlying connectivity pattern 
affects the dynamical interactions among the examined 
elements [TJ [2l [3] . Going beyond characterizing the net- 
work topology with the essential degree distribution [4], 
extensive research has focused on the influence of degree- 
degree correlations [2 . The degree of a node equals to 
the total number of links attached to it. A network with 
a positive degree-degree correlation is usually called as- 
sortative, and implies that nodes with a similarly small or 
large degree tend to be connected to each other [21 [5] . If, 
by contrast, the nodes tend to be connected to nodes with 
a considerably different degree, the network is called dis- 
assortative, referring to a negative degree-degree correla- 
tion. Assortativity is typically found in social networks, 
disassortativity in biological and technical networks [6]. 
The impact of correlations on the dynamics appears to 
be non-trivial [7 and has so far been discussed by analyz- 
ing specific processes only. In this respect, assortativity 
seems to be a suppressing factor for the spatio-temporal 
spreading of epidemic diseases [8 , while at the same time 
disassortativity has been suggested to prevent the propa- 
gation of perturbations in protein networks [9] and to en- 
hance the robustness of declining company networks [10] . 

In this paper, we show that many spreading pro- 
cesses can be categorized into two fundamentally differ- 
ent regimes, in which either assortativity or disassorta- 
tivity has a decelerating effect. For this aim, we broadly 
define a spreading process as the propagation of a certain 
nodal state along the links of a network. The two regimes 
are then given by the probability of a node to respond 
to the current states of its neighboring (i.e. connected) 
nodes. In regime I the probability of adopting the prop- 
agating state is lower for nodes with a large degree than 
for nodes with a small degree, given that the same ratio 



of neighboring nodes already is in that state (see illustra- 
tion in Table |l|. In regime II the probability is higher for 
nodes with a large degree. Such a connectivity dependent 
local response mechanism is found in existing models for 
various phenomena, even though they may strongly differ 
besides in their specific features. An example of regime I 
is a model for a declining company network, whereas the 
probability for a company disappearing is inversely pro- 
portional to its degree Hence, after loosing the same 
ratio of connected firms a company with a large (initial) 
degree has still more connections and thus a lower prob- 
ability to disappear than a company with a small degree. 
An example of regime II is the spreading of epidemics 
on the worldwide air-transportation network [TT], since 
a highly connected city is more susceptible than a city 
with less connections, given that the same ratio of adja- 
cent cities has an equally infected population. Further 
examples of categorizable models are listed in Table |lj 

TABLE I: (Color online) Examples of models for dynamic pro- 
cesses and their categorization based on the local response of 
a node to the states of its neighboring nodes. In the illustra- 
tion node A changes its state with probability Pa and node B 
with Pb respectively, whereas the same ratio of neighboring 
nodes is already in the propagating state (red (dark) color). 



Regime /: Pa < Pb 
Declining company network [10] 
Extinction of species [12 ^ 



Regime II: Pa > Pb 
Reaction-diffusion processes [13] ^ 
Global epidemics 11 
Dissemination of information [14] ^ 
Cascading failures in power grids [T5] 



^Similarly to the declining company network (see text) an ex- 
tinction sequence follows the deletions from the least to the most 
connected species (being a deterministic model in fact). 

node with a large degree has a higher probability to receive an 
active particle than a node with a small degree, given that the same 
ratio of neighboring nodes has the same density of active particles. 

^Large-degree nodes are likely to get the information at a lower 
ratio of neighboring nodes being already informed. 

^If a certain ratio of neighboring nodes fails, a node with a large 
degree has a higher probability to become overloaded. 
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In order to focus on these two local response mech- 
anisms and to analyze the regime dependent effect of 
degree-degree correlations on the resulting spreading dy- 
namics, we draw upon the illustration of a decaying undi- 
rected network. The spreading is captured by a minimal- 
istic two-state model, whereby the nodes are either fully 
functional or failed. Let us stipulate that a node becomes 
stressed as one of its neighboring nodes fails and that the 
stress, in turn, increases the probability that this node 
fails during the next time step [16 . It is then natural to 
assume that the stress Si on a node i increases with the 
ratio Xi = k( /ki of failed neighboring nodes k( to the 
initial node degree, ki. We further assume that 5^ = if 
k{ = and that the maximum stress 5^ = 1 is indepen- 
dent of the inital degree and is reached when k{ = k^. 
In order to vary the relative influence of the initial node 
degree we introduce the following formula for the stress 
Si{xi): 



Xikj 



(1) 



where 6 is the response parameter. As shown in Fig. [l] 
the stress equals to if 6> = 0. For 6> < we obtain 
regime I, as the absolute number of functional neighbor- 
ing nodes becomes more important. For example, having 
lost 50% of the neighbors, a node i initially with ki = 30 
neighbors is less stressed than a node j with an initial 
degree of kj = 10. Analogously, 6> > corresponds to 
regime II, as the absolute number of failed nearest neigh- 
bors becomes more important. For ^ the loss of just 
one neighbor causes maximum stress, independently of 
ki. In the limit ^ <C the stress is 5^ ~ if < < 1, 
but jumps to = 1 when all neighboring nodes have 
failed. It is worth remarking that Eq. ([T]) can be re- 
placed by other functions which qualitatively reproduce 
the two regimes. 
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FIG. 1: (Color online) Illustration of the stress function de- 
scribed by Eq. [l] For ^ < the response of nodes to the 
loss of their neighboring nodes is according to regime I and 
inversely, for ^ > 0, we obtain response regime II. 

Starting at time t = with Nq nodes, the network 
decay is governed by both random (i.e. independent) 



and stress induced node failures, formulated by the time 
dependent failure rate of a node i: 



= p^qsi{xi(t)), 



(2) 



where Fiit) — fi{u)du is the life time distribution 
function. The parameters p > and q > account for 
the random failures and for the influence of the stress re- 
spectively. The relation to different classes of spreading 
models is given by the parameter range of Eq. ([T]) and 
Eq. ([2|. For instance, putting p ^ in the limit ^ ^ 
we readily recover the susceptible-infected (SI) model as 
applied to study epidemics in which individuals are either 
in the susceptible or infected state [17] . The proportional 
increase of A^, as ^ = 0, corresponds to the Bass diffu- 
sion model [18 . A further example are binary threshold 
models as applied for social contagion processes |19l |20] ; 
assuming a high (low) threshold these models can be ap- 
proximated by setting a correspondingly low (high) value 
for ^ at p ^ and a high value of q. 

We studied the spreading processes on finite size scale- 
free networks, differing in their global level of correlation. 
The degree distribution of scale-free networks follows a 
power law, P{k) ex which is found in many real life 
systems, whereas the characteristic degree exponent usu- 
ally lies in the range 2 < 7 < 3 [1 . The global level 
of correlation is commonly quantified by the Pearson co- 
efficient r [21 . This measure is defined in the interval 
[—1,1] whereas r = corresponds to an uncorrelated net- 
work and a positive (r+) (negative (r_)) value denotes 
positive (negative) correlation. We first constructed un- 
correlated networks with a given scale-free exponent 7 
according to the algorithm presented in [22 , restricting 
the degree ki of each node i to m < ki < \^No with 
ki being even. We subsequently applied the reshuf- 
fling method [23] to produce networks with the desired 
correlation coefficient, being in the range —0.3 < r < 0.5. 
The positive range includes higher absolute values than 
the negative range in order to better account for real 
life networks. Their maximum degree is not restricted 
to ViVo, so that a higher positive correlation below this 
structural cut-off might become concealed in the value of 
r as nodes with a larger degree naturally induce negative 
correlations [22] . 

The velocity of the stress propagation determines the 
nodal life expectancy, which represents the average time 
span until a randomly chosen node fails: the slower the 
spreading, the larger its value. Estimating this measure 
by extensive Monte Carlo simulations shows evidence 
for disassortativity decelerating spreading processes of 
regime I and assortativity decelerating those of regime 
II. This result is presented in Fig. [2| where we com- 
pare the nodal life expectancies T(r+) in assortative and 
r(r_) in disassortative networks, being normalized with 
the corresponding value of the uncorrelated case, r(0). 
The scale-free networks have Nq = 10^ nodes, character- 
istic exponent 7 = 2.5 and minimum degree m = 2. The 
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parameters controlling the nodal failure rate are set to 
p 0.0001 and (7 = 1 - j9 = 0.9999. 



1.6 
1.4 
1.2 



0.8 
0.6 
0.4 
0.2 



regime I 
oo 



regime II 

^AAAAa^Aa^^^^^^AAaI^ 



AO 



g 



r = 


a J 


H-B r = 


0.4 


A A r = 


0.5 


A A r = 


-0.1 


□ □ r = 


-0.2 




-0.3 



'S|n^ „fA 



e 

FIG. 2: (Color online) Relative nodal life expectancy 
T{r)/r{0) in scale-free networks with A^o=10^, 7 = 2.5 and 
m = 2 for different values of r versus the response parameter 
0. The inset shows the absolute values for t(0) in the uncor- 
related reference networks. The dotted lines serve as a guide 
to the eye. 

Starting in regime I with a very low susceptibility to 
failures, i.e. ^ <C 0, the decay is dominated by random 
failures as the stress propagation is slow (see inset to 
Fig. [2|. Increasing strongly increases the influence of 
Si on the nodal life expectancy, with r(r+) < r(0) < 
r(r_). These characteristics are more pronounced for 
larger values of |r|. After observing T(r+) ~ r(0) ~ 
r(r_) in the intermediate range, the spreading sustains 
with high absolute velocity in regime II at ^ > 0. Now 
the nodal life expectancy in assortative networks becomes 
larger, with r(r_) ~ r(0) < T(r+). Finally, in the limit 
^ ^ the decelerating effect of the assortativity seems to 
persist, being consistent with earlier findings on epidemic 
spreading [8]. Although we are focusing on scale- free 
networks we may expect qualitatively similar results for 
other correlated networks. 

The observed behavior can be attributed to the role of 
the nodes with a large degree and their distribution in the 
network. In regime I they retain the propagation of the 
stress (see Fig. [T]), acting as delayers^ whereas in regime 
II they become accelerators. In assortative networks they 
become clustered in regions of high connectivity, which 
can be identified for example by the /c-core decomposi- 
tion [24l [25] . A k-core is the maximum subgraph with all 
nodes having minimum degree /c, and a /c-shell contains 
the fraction of nodes belonging to the /c-core but not to 
the (/c + l)-core, see Fig. [3](a) for illustration. Hence, in 
case of regime I the delayers become bunched into higher- 
order /c-shells, allowing for a fast stress propagation in 
the lower-order /c-shells (see Fig. [3](b)). Analogously, in 
the case of regime II the accelerators become separated, 
effectively decelerating the spreading velocity within the 
lower-order /c-shells (see Fig. |3](c)). In this regime, the 
general spreading progression from nodes with a large 
degree towards nodes with a small degree is again con- 



sistent with epidemic models [26]. Due to the degree 
distribution of the scale- free networks the vast majority 
of the nodes remains in lower-order /c-shells (see inset to 
Fig. [3] (c)) so that this opposite effect becomes directly 
reflected in the average nodal life expectancy. 
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FIG. 3: (Color online) (a) Illustration of the /c-core decompo- 
sition of a network before (left) and after (right) reshuffling 
towards a higher positive correlation. Each shaded area corre- 
sponds to a different /c-shell. (b)(c) Life expectancy {rk- shell) 
of the /c-shells for assortative (r = 0.5, inversed triangles) and 
disassortative (r = —0.3, squares) scale- free networks with 
ATo^lO^, 7 2.5 and m = 2. (b) In regime I with 6* = -2 
and (c) in regime II with — 2. All other spreading model 
parameters are as of Fig. [2] The inset illustrates the relative 
size of the /c-shells. 

In regime I the retention times of the stress in the de- 
layers are long, so that the propagation, when possible, 
bypasses through nodes with a smaller degree. The op- 
posite is true for the accelerators in regime II. Hence, the 
different encapsulation of the delayers (respectively ac- 
celerators) as captured by the /c-core decomposition also 
changes the overall spreading paths through the entire 
network and thus the time needed until the stress hits 
a randomly chosen node. To illuminate this effect, we 
examined the 'efficient paths' through which the stress 
propagates most likely. Its basic concept has been in- 
troduced in [27 for optimal routing in communication 
networks. We first calculated the length of a path Vi^j^ 
connecting node i with node j and containing the set of 
nodes V-p, as 



(3) 



where is a parameter accounting for the degree depen- 
dent retention times. The efficient path length is then 
given by the minimum value of L'^{Vij;iy) for all pos- 
sible paths between nodes i and j. Averaging over the 



4 



efficient path lengths between ah possible pairs of nodes 
eventually leads to the average efficient path length (/^). 
If = 0, we retrieve the geodesic shortest path. The re- 
sulting average efficient path lengths for different values 
of V are depicted in Fig. |4] 
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FIG. 4: (Color online) Average efficient path length of cor- 
related networks, (/^(r)), normalized with the uncorrelated 
reference case, (/^(O)), versus the parameter u. Each curve 
corresponds to an average over 100 independent realizations 
of networks, having parameters as of Fig. |2] 

The values for (P) indeed show a qualitatively sim- 
ilar behavior as the result of the stress propagation 
(see Fig. [2|. For u < disassortative networks ex- 



hibit a larger value for (/^). The differences at ~ 
stem from the fact that assort at ive networks have a sig- 
nificantly larger (and disassortative networks a slightly 
larger) geodesic shortest path length than the uncorre- 
lated reference networks. At > the average effi- 
cient path length becomes larger for assort at ive networks. 
Given this good agreement, the stress propagation seems 
to follow the efficient paths, suggesting that (/^) might 
be a robust indicator for the relative influence of degree- 
degree correlations on the spreading velocity. It further 
supports our hypothesis, which is to categorize spreading 
processes into the two proposed regimes. 

In summary, we have drawn a global picture on how 
degree-degree correlation affects the velocity of spreading 
processes. If the propagation follows regime I, disassor- 
tativity decelerates the velocity. Conversely, in regime II, 
it is assortativity which prevents a fast spreading. This 
finding and the underlying role of the nodes with a high 
connectivity are believed to be important for either im- 
proving the robustness of networks against undesirable 
spreading processes such as epidemic diseases or to ac- 
celerate the desired propagation of, for example, infor- 
mation. 
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